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ABSTRACT 

The gas-phase chemistry of water in protoplanetary disks is analyzed with a model 
based on X-ray heating and ionization of the disk atmosphere. Several uncertain pro- 
cesses appear to play critical roles in generating the column densities of warm water that 
are detected from disks at infrared wavelengths. The dominant factors are the reactions 
that form molecular hydrogen, including formation on warm grains, and the ionization 
and heating of the atmosphere. All of these can work together to produce a region of 
high water abundances in the molecular transition layer of the inner disk atmosphere, 
where atoms are transformed into molecules, the temperature drops from thousands to 
hundreds of Kelvins, and the ionization begins to be dominated by the heavy elements. 
Grain formation of molecular hydrogen and mechanical heating of the atmosphere can 
play important roles in this region and directly affect the amount of warm water in 
protoplanetary disk atmospheres. Thus it may be possible to account for the exist- 
ing measurements of water emission from Tauri disks without invoking transport of 
water from cooler to warmer regions. The hydroxyl radical OH is under-abundant in 
this model of disk atmospheres and requires consideration of additional production and 
excitation processes. 

Subject headings: stars:planetary systems: protoplanetary disks - X-rays:stars - astro- 
chemistry 



1. Introduction 

The water molecule plays a significant role in the evolution of protoplanetary disks and in 
the formation of small and large bodies including planets (e.g., Ciesla & Cuzzi 2006; Jewitt et 
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al. 2006). Recent results from the Spitzer Space Telescope (Spitzer) reveal rich mid infrared (MIR) 
water emission spectra from 10-35 microns that appear to be a common characteristic of T Tauri 
disks (Carr & Najita 2008; Salyk et al. 2008). Assuming thermal equilibrium populations, typical 
temperatures and column densities have been estimated to be T ~ 500 — 1000 K and 7V(H20) ~ 
10 17 — 10 18 cm~ 2 , with emitting areas of the order of a few AU in radius. Mid infrared emission 
from other molecules (HCN, C2H2, OH, CO2) is also detected with similar temperatures and likely 
similar emitting radii, e.g., Carr and Najita (2008). These MIR lines appear to probe lower gas 
temperatures and larger disk radii than the previously known suite of gaseous inner disk diagnostics. 
UV fluorescent emission from H2 and near-infrared transitions of CO probe warm gas temperatures 
(1000-3000 K) and primarily the inner < 1 AU of the disk (Najita et al. 2007). Measurements of 
molecular line emission have been obtained for many other protoplanetary disks by Spitzer and 
should become available for analysis and interpretation in the near future. 

Our earlier model of the thermal-chemical structure of an X-ray irradiated generic T Tauri 
disk atmosphere (Glassgold et al. 2004; henceforth GNI04) included the synthesis of H2O by 
temperature-sensitive neutral radical reactions (Table 1 of that paper), but it did not predict 
strong water emission. Possible explanations for this deficiency might be that water is dredged up 
from deeper parts of the disk atmosphere by turbulent mixing, or that it is transported radially 
by inward migration of intermediate-sized bodies and desorbed (Ciesla &; Cuzzi 2006; Najita et 
al. 2007; Salyk et al. 2008). In this case, measurements of water abundances in disk atmospheres 
might be used to trace dynamical processes such as turbulence and migration. However, disk gas 
phase chemistry is affected by several often poorly known physical processes, so it is important to 
fully explore non-dynamical explanations. In this paper, we analyze some of these processes in light 
of the recent Spitzer H2O observations by exploring alternative or modified chemical pathways that 
may enhance the abundance of water and other molecules in disk atmospheres. A longer term goal 
is to attempt to use observations of water to infer the physical properties of protoplanetary disks. 

This study complements recent work on disk chemistry by focusing on the thermal-chemical 
properties of an X-ray irradiated disk atmosphere that has experienced grain growth and settling. 
X-ray induced chemistry has a long history going back 25 years (Krolik and Kallman 1983). Ref- 
erences to early work is given in the pioneering papers by Lepp and Dalgarno (1996) and Maloney 
et al. (1996). The latter calculate the water abundance in X-ray-transition regions, as do Mei- 
jerink and Spaans (2005). Stauber et al. (2005) develop X-ray chemistry for spherically symmetric 
envelopes around YSOs with specific emphasis on water (Stauber et al. 2006). 

Many disk chemical models consider the outer region of disks, whereas we concentrate on the 
inner < 10—20 AU. Several recent papers treat inner disk gas phase chemistry (Markwick et al. 2002; 
Nomura and Millar 2005; Nomura et al. 2007, 2009; Agundez et al. 2008; Gorti and Hollenbach 
2008; Woods and Willacy 2009). Nomura et al. (2007) focus on H 2 , and Nomura et al. (2009) follow 
molecules along prescribed stream lines and find that the abundances of many organic molecules 
change as they enter the inner disk region, notably methanol but not water. Markwick et al. and 
Agundez et al. include X-ray ionization in terms of a phenomenological ionization rate, whereas 
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Woods and Willacy use a more complete theory based on Gorti and Hollenbach (2004). Agundez 
et al. carry out time-dependent abundance calculations, in some cases at fixed temperatures. We 
shall discuss the water abundances obtained by these authors in the Discussion section below. 

In this paper we analyze a number of physical processes that can affect the abundance of water 
in the warm part of the disk atmosphere (T > 200 — 300 K) using the GNI04 model. In Sec. 2, we 
discuss the role played by dust through the thermal coupling of the gas and dust; the efficiency of 
the gas-phase reactions that produce H2 and H2O; the role of H2 and H2O formation on dust grains; 
and finally the effects of enhanced mechanical heating of the atmosphere. Implementation of these 
processes with our thermal-chemical model, along with new information about their efficiency, 
indicates that column densities of warm water in the range observed can be obtained without 
invoking radial or vertical transport. The implications of this results are taken up in Sec. 3 and in 
the Conclusions in Sec. 4. 

2. Thermal-Chemical Model of X-ray Irradiated Disks 

2.1. The Basic Model 

This work builds on our earlier model of the thermal-chemical structure of the inner disk of 
a T Tauri star (GNI04; see Meijerink et al. 2008 for corrections and additional applications). The 
disk is irradiated by stellar X-rays which ionize and heat the disk surface. The X-ray luminosity 
is Lx = 2 x 10 30 ergs _1 ; a thermal spectrum is used with Tx = 1.0 keV and a low-energy cutoff of 
Eq = 200 eV. As in GNI04 and our earlier work, cosmic rays are ignored on the basis of their being 
blown away from the inner disk by the stellar wind. The total gas density distribution and the 
dust temperature are taken from D'Alessio et al. (1999, 2001), and the gas-phase temperature and 
abundances are solved for self-consistently keeping the density and dust temperature fixed. In a 
more general model that treats the dust and gas separately, a new equilibrium would be calculated, 
including the effects of changing the gas temperature on the hydrostatic balance and on the dust 
temperature. 

The gas temperature in the current model is determined by balancing heating from X-rays and 
mechanical processes with cooling by line radiation. GNI04 employed a simple chemistry with 25 
species and 125 reactions that included the heavy atoms, C and O, as well as H and He. In order 
to provide a more complete theory of X-ray ionization, the chemical code used here treats all of 
the most abundant heavy atoms and includes nearly 400 reactions (Meijerink & Glassgold 2009). 
As mentioned in the Introduction, chemical synthesis is achieved by neutral radical reactions with 
moderate barriers that can be overcome in the warm upper atmosphere of the disk. Ionic reactions 
also contribute, but their most important role is to destroy molecules. One difference with GNI04 is 
that we use depleted abundances, defined relative to hydrogen (nuclei) as usual, with xq, = 3.5 x 10~ 4 
and xc = 1.4 x 10~ 4 . The thermal part of the program is unchanged from GNI04 even as changes 
are made in the chemistry. For example, Nell and water cooling are still ignored, although their 
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inclusion could change the quantitative aspects of the present results. In addition, a significant 
fraction of the stellar NIR flux may be captured by water molecules. A consistent theory of the 
thermal effects of water in disk atmospheres requires consideration of both heating and cooling. We 
postpone treating this important topic to a future study. Also omitted in this pure X-ray irradiated 
disk model are the effects of UV radiation, which can affect the abundances of OH, H2O and other 
species. Efforts are under way to include both UV radiation and X-rays in our models. 

Dust grains play a direct role in the thermal balance via the familiar collisional coupling of the 
dust and the gas. Cooling of the gas by the dust becomes important in the region we refer to as 
the transition region where the temperature and ionization level are rapidly decreasing and atomic 
species are transformed into molecules. Since the chemistry is sensitive to temperature, the dust 
indirectly affects the chemical abundances. In the GNI04 model, it was assumed that the dust in 
the inner disk had experienced growth and settling. The grain area per H nucleus was given by 
(Eq. A4 of that paper), 

n d ^ 2 . Pd/Pg rn 1/2 

— < a >= r: , flgcom = {aia 2 ) 1 , (1) 

m a g eom (4vr/3)p 

where ua and pa are the volumetric number and mass density of the dust grains, p g is the mass 
density of gas, p is the internal density of a typical grain, m is the mean mass of the gas per H 
nucleus, and a\ and ai are the minimum and maximum grain radii for a MRN distribution. The 
key parameter in Eq.[T]is the dust to gas ratio Pa/ Pg divided by a ge0 m> the geometric mean of a± and 
a,2- In GNI04, these parameters were set at Pa/ Pg = 0.01 and a ge0 m = 0.707 pm, following Wood et 
al. (2001). This choice corresponds to a factor of 20 increase in QWom and a factor of 20 decrease in 
effective grain surface area relative to the interstellar medium. Wood et al. actually used a MRN 
distribution with a\ = 0.01 pm, an exponential cutoff at a = 50 pm and 02 = 1 mm. Had we used 
the grain model in D'Alessio et al. (2001; a% = 0.005 pm and 02 = 1mm), then a ge0 m would be 
2.24 pm, or 3.2 times larger than in GNI04. In this paper we follow GNI04 but also consider larger 
values of a geom , i.e., a geom = 7.07 pm as well as a geom = 0.707 pm. This range gives grain areas per 
H nucleus that are roughly consistent with the dust depletion deduced by Furlan et al. (2006) for 
the upper atmospheres of T Tauri disks observed with Spitzer. 



2.2. The Dependence of the Water Abundance on the Physical Properties 

Only gas phase reactions were considered by GNI04 for the formation of H2 because the grains 
in disk atmospheres were considered too warm for the usual grain synthesis familiar from the 
interstellar medium (e.g., Hollenbach & Salpeter 1971) to be efficient. Instead, H2 was formed by 
radiative attachment of electrons 

e + R^R' + hv H~ + H^H 2 + e, (2) 

and by three-body reactions, 

3H^H 2 + H 2H + H 2 ^2H 2 , (3) 
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the latter accompanied by collisional dissociation, 

H + H 2 ^3H H 2 + H 2 ^H 2 + 2H. (4) 

These gas phase reactions led to H 2 abundances approaching ~ 10~ 3 in the surface region, a level 
sufficient to produce significant amounts of warm CO (~ 10 17 — f 18 cm" 2 ) by a sequence of neutral 
radical reactions that begins with the slightly endothermic reaction 

O + H 2 ^H + OH (5) 

and leads to water by the exothermic reaction, 

OH + H 2 ^H 2 + H. (6) 

Both of these reactions have barriers that require warm conditions to be effective, i.e., T > 200 — 
300 K. OH also leads to 2 and CO by fast exothermic radical reactions without significant barriers, 

OH + 0^0 2 + H C + OH^CO + H. (7) 

OH and H 2 are mainly destroyed by charge transfer with H + , 

H+ + OH^OH+ + H, (8) 

H + + H 2 O^H 2 0+ + H. (9) 

The chemistry defined by the above equations starting with Eq. 5 can be approximated by the 
following formula for the abundance of water, 

x(H 2 0) / x(H 2 



<n 9^) do) 



z(O) v ; v^( H+ ) 

where 

K = „ 10 -18 T 4.27 exp( _ 482 3/T). (11) 

kgkg 

The ratio of water to atomic oxygen depends quadratically on the ratio of H 2 to H + because 
the synthesis of water proceeds via two radical reactions with H 2 , (Eqs. 5 and 6), and because 
the precursor radical OH as well as H 2 are destroyed primarily by fast charge exchange with 
H + (Eqs. 8 and 9). The factor n(T) is sensitive to temperature because the product of the rate 
coefficients, k$ and k$, for the two synthesis reactions is temperature-sensitive, whereas the product 
of the two rate coefficients, ks and kg, for the two ionic destruction reactions is not. Equation 10 
ignores several factors that can limit its accuracy: the role of the backward reactions to Eqs. 5 and 
6 at high temperatures; the contribution of ion-molecule reactions to the synthesis of water; and 
the destruction of OH and H 2 by He + ions. 

The main purpose of introducing Eq. 10 is to elucidate the underlying physics of the molecular 
transition and not to reproduce the full model calculations. The main conclusion from this ap- 
proximate treatment is that the water abundance in disk atmospheres is sensitive to two quantities 
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that vary rapidly in the transition zone between the very warm upper atmosphere and the cool 
near-midplane region, the temperature and the ratio of the molecular hydrogen abundance to the 
H + abundance. Any process that affects these quantities can change the abundance of water. For 
our model of an X-ray irradiated disk, the X-rays play a key role since they heat and ionize the 
gas. 

The GNI04 model calculations yielded small column densities of warm water (< 10 13 cm" 2 at 
T > 200 — 300 K). In contrast, observations indicate the presence of abundant warm water, e.g., 
Carr et al. (2008) report iV(H 2 0) = 6.5 x 10 17 cm" 2 at T = 500 - 600 K in AA Tau, and Salyk 
et al. find similar columns in AS 205 A and DR Tau. According to our calculations, the water is 
probably located in the vertical thermal-chemical transition region that occurs between the hot and 
cool parts of the disk, where the gas temperature is relatively high, e.g., in the 300-2000 K range. 
The observed water lines are in emission, which further suggests that the lines form above the MIR 
dust photosphere. The deficiency of the GNI04 model in this respect suggests three possibilities: 
(1) other processes heat the atmosphere over larger column densities, so that the water deeper 
down in the atmosphere is put into emission; (2) water is transported dynamically and mixed into 
the upper atmosphere where it is heated; (3) more efficient chemical pathways form water in the 
disk atmosphere. We mainly focus on options (1) and (3) in this paper. 



2.3. Gas Phase Reactions 

Reaction [5] is an important destruction mechanism for H2. The recent critical review by 
Baulch et al. (2005) contains laboratory measurements of this reaction from 300-3500 K. The low- 
temperature data (T > 300 K) are 20 years old and were obtained with shock tubes or flash 
photolysis. GNI04 used the value recommended in the earlier Baulch et al. (1992) review, 

k 5 = 8.5 x 10- 20 T 2 - 67 e-s^cmV 1 . (12) 

The more recent experiments all refer to the higher-temperature regime, T = 1500 — 3500 K, and 
the above fit is still satisfactory at lower temperatures. Nonetheless, considerable uncertainty must 
be attached to Eq. [5] (~ 50 — 100%), and even more when it is extrapolated below 300 K. GNI04 
mistakenly used the exponent 1.67 instead of 2.67, seriously underestimating the production of 
oxygen molecules while underestimating the destruction of H2 via Eq. This error was corrected 
in MGN08, but even with this change only a small amount of warm water was generated, e.g., 
iV(H20) < 10 13 cm" 2 at 1 AU. Significant water abundances, say x(H20) > 10" 6 , were only 
achieved at large vertical column densities (Nr > 10 22 cm" 2 ) where the dust temperature is barely 
above freeze-out. This situation is illustrated by the dotted curves in the middle and bottom panels 
of Figure 1 for a radial distance of 1 AU. 

Another critical gas phase process is three-body formation of H2, the reactions in Eq. [3j In the 
past, these three-body rate coefficients were obtained with detailed balance from measurements of 
the inverse reactions, usually those of Jacobs et al. (1967). Flower and Harris (2007) have used new 
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equilibrium constants to obtain three-body formation rate coefficients that are significantly larger 
than those given in the review of Cohen and Westberg (1983) and used by GNI04. For example, 
the first forward reaction in Eq. 3 is, as given by Flower and Harris, 

fc 3 (H^H 2 + H) = 1.44 x io- 26 T" L54 cmV 1 , (13) 

which is to be compared with the Cohen and Westberg rate, 8 x 10~ 33 cm 6 s _1 . 

The top panel of Figure 1 shows the change in the abundances of hydrogenic species obtained 
using the new rate coefficients. The difference between old and new rate coefficients is not manifest 
until Nu > 5 x 10 20 cm -2 because the H2 formation proceeds mainly via the H~ route, Eq. [21 at 
high altitudes. The strong rise in the H2 abundance as A/h approaches 5 x 10 20 cm -2 occurs because 
the three-body rate increases rapidly with decreasing temperature in a region where the atomic H 
density is also increasing towards 10 10 cm -3 . As a result, the H-H2 transition to be is shifted to 
smaller vertical column densities, which means to warmer regions higher in the atmosphere. As can 
be seen Figure 1, this also leads to a significant increase in the amount of warm water, aided by a 
drop in the abundance of H + according to Eqs.[8]and[H The incomplete transitions of atomic O to 
OH and H 2 (to abundances ~ 1CT 6 ) now occur in the same general region as the hydrogen and 
carbon transitions. The solid curves in Figure 1 define the first in a series of comparison models to 
be introduced in this section ("model 1"). The other models are defined in Table 1; their specific 
properties are given later in Table 2. Note that the parameter ag eom listed in Table 1 determines 
the dust cooling of the gas as well as the formation of H 2 on grains. 



2.4. The Formation of H 2 on Dust Grains 

As discussed in Sec. 2.2, we previously considered only gas phase pathways for forming H 2 on 
the assumption that dust grains in disk atmospheres are too warm for efficient grain synthesis. This 
idea is well founded in the standard theory of H 2 formation in interstellar clouds (e.g., Hollenbach 
and Salpeter 1971) and substantiated by low-temperature laboratory experiments using interstellar 



Table 1:: Chemical Models 



Model 


Chemistry 


ageom [/^m] 


a h 


1 


pure gas phase 


0.707 


0.01 


2 


H 2 formation on grains 


0.707 


0.01 


3 


H 2 formation on grains 


7.070 


0.01 


4 


H 2 & H 2 formation on grains 


0.707 


0.01 


5 


pure gas phase 


0.707 


1.00 


6 


H 2 formation on grains 


0.707 


1.00 
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Fig. 1. — : Hydrogen (top) and oxygen (middle) species at 1 AU plotted vs. vertical column 
density for old (dotted line) and new (solid line) three-body reaction rate coefficients discussed in 

GNI04 (dotted line) and this paper (solid line). See Sec. 2.3 for details. With the new reaction 
rates, the H-H2 transition, where x(H) = x(H2), is shifted upwards by a factor of three in vertical 
column from JVh = 7 x 10 21 cm -2 to Nr = 2 x 10 21 cm -2 . The transition to water is also similarly 
shifted upwards. The solid curves represent model 1 in Tables 1 and 2. 
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dust analogs and interpreted with theoretical models based on the physical adsorption of H atoms 
(e.g., Vidal et al. 2004; Perets et al. 2005, 2007). The efficiency of H2 formation is close to unity 
in a narrow range (width ~ 5K) of temperature near 15 K, depending on the experimental sample. 
At lower temperatures, incident H atoms are insufficently mobile, and at higher temperatures 
they leave the surface before forming molecules. Invoking a model based on chemical adsorption, 
Cazaux and Tielens find that H2 can form with moderate efficiency (< 0.2) on warm dust with 
temperatures up to 900 K (Cazaux & Tielens 2004; Cazaux et al. 2005). We adopt this model 
here because it is a reasonable extension of existing understanding and because it supported by 
laboratory experiments with graphite at intermediate temperature of a few hundred Kelvins (e.g., 
Zecho et al. 2004; Hornakaer et al. 2006a, b). 

We add grain formation to the program described earlier by writing the H2 formation rate per 
unit volume in the standard form, i.e., 1/2 times the destruction rate of atomic H: 

R= I n (Il>(H)n d < vra 2 > eS(T g ,T d ). (14) 

The two factors after the 1/2 give the incident H atom flux, the next two specify the grain surface 
area per unit volume, and the last two are the product of the formation efficiency and the sticking 
probability (or, more simply, the overall efficiency). We use GNI04 for the mean grain area per 
unit volume, based on the MRN distribution and an internal dust grain density of 3 gcm~ 3 , to 
obtain (in units of cm 3 s~ 1 ), 

R = 4.32 x lO- 19 ^ n( H)n H eS x (^) (15) 

We use the right panel of Figure 2 of Cazaux et al. (2005) for the formation efficiency e as a function 
of the dust temperature Td for the one case they considered where the gas and dust temperatures 
are unequal, T g = Trf + 500K. A similar gaseous temperature enhancement is expected in the region 
in which the molecular emission is likely to arise. We estimate S = 0.1 from Burke & Hollenbach 
(1983). 

The upper panel of Figure 2 gives the results for the hydrogen family at R = 1 AU. Three sets 
of abundance curves are plotted vs. vertical column density for three levels of grain formation: no 
grain formation (dotted line); grain formation with a geom = 7.07;Um (dashed curve); grain formation 
with a g eom = 0.707)Um (solid curve); they correspond respectively to models 1, 3 and 2 in Table 1. 
Note that the usual interstellar value is a g = 0.035 ^m. In accord with the discussion in Sec. 2.1, the 
two choices of a ge0 m correspond to reduction factors of 20 and 200 in the grain area per H nucleus 
compared to ISM values. If we define the location of the H to H2 transition by the vertical column 
where x(H) = x(H2), then the transition occurs higher in the atmosphere (at smaller vertical 
columns) as the grain area increases. For the case a ge0 m = 0.707//m, the H-H2 transition, where 
x(H) = x(H2), is shifted upwards by a factor of four in vertical column from Ah = 2.4 x 10 21 cm -2 
to Ah = 6 x 10 20 cm -2 . The bottom panel of Figure 2 plots the temperature for the three cases in 
the upper panel. The dotted and solid temperature curves are roughly the same because they have 
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the same grain area, i.e., a ge0 m = 0.707/im. For the dashed curve, a ge0 m = 7.07^m, the temperature 
is raised above the dust temperature because the gas and dust become less-well coupled thermally 
as the grain surface area decreases. Comparing the upper with the lower panel of Figure 2, we see 
that the H-H2 transition for the solid curves (a ge0 m = 0.707/im) occurs roughly in the middle of 
the transition region and that this model yields a significant column of warm water (T > 300 K), 
N(R 2 0) «2x 10 17 cm" 2 . 

Figure 2 shows that grain formation of H2 is important in the warm transition region of the 
disk atmosphere. Increasing the abundance of H2 also increases the abundances of OH and H2O, 
which leads to a decrease in the abundance of H + through the reactions in Eqs. [8] and [9j This 
adds a degree of runaway to the abundances since decreasing H + increases the amounts of OH 
and H2O. The two main oxygen species, O and H2O, are plotted in the middle panel of Figure 
2. The increased abundance of H2 manifests itself in a significant enhancement of H2O in the 
warm transition region. Without grain formation of H2, full conversion of O to H2O does not 
occur until a vertical column Ah = 1-5 x 10 23 cm -2 has been reached, at which depth dust cooling 
has overpowered gas heating and the gas and dust temperatures are the same. At R = 1 AU, this 
temperature is of order 100 K, so such a model cannot account for the warm water observed at near- 
and MIR wavelengths. However, with grain formation of H2, full conversion occurs higher in the 
atmosphere where the gas temperature is elevated above that of the dust. For example, for the case 
Qgeom = 7.07/im, the O-H2O transition occurs near ATjj = 10 cm and, for Ogeom = 

0.707/iin, 

it occurs at A% = 4 x 10 20 cm -2 . In both cases, there is a warm layer where almost of all 
the available oxygen is in water (recall that x<j = 3.5 x 10 -4 and xq = 1.4 x 10~ 4 ). For the 
case Ogeom = 0.707/xm, this layer contains a water column, A^^O) ~ 2 x 10 17 cm -2 , and for 
Ogeom = 7.07/xm, the water column is A^^O) « 4 x 10 17 cm~ 2 . This increase occurs because 
the increase in a ge om implies a decrease in dust cooling, which results in a shift downwards of the 
warm region into denser gas. These columns of warm water approach the values obtained from 
recent Spitzer observations assuming thermal populations, > 10 17 cm -2 (Carr & Najita 2008; Salyk 
et al. 2008). The OH abundance also rises as the transition to water is approached, reaching 
abundance levels ~ 10 -7 — 10 -6 and then decreasing rapidly with increasing vertical column of 
hydrogen. 

Below these layers of near maximum water abundance high in the atmosphere, the water 
abundance drops two orders of magnitude before recovering at very large depths. The origin of this 
abundance swing can be found in the approximate formula Eq. [TO) Once H2 has been fully formed, 
the H2O/O ratio is governed by the relative importance of the temperature decline, as the gas and 
dust temperatures become more and more closely coupled, and by the decrease in the abundances 
of the ions H + and He + which destroy H2O. With increasing depth, the decrease in temperature 
first reduces the production of water, but then the attenuated X-ray flux (and consequent low ion 
abundance) removes any effective destruction so that its abundance increases again. 

The full conversion to water is closely associated with the two other chemical transitions, H 
to H2 and C to CO. These transitions are displayed in Figure 3 at two smaller radii, R = 0.25 and 
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Vertical column density [cm 2 ] 



Fig. 2. — : Main hydrogenic species, atomic oxygen and water, and gas temperature at 1 AU 
plotted vs. vertical column density for three cases of grain formation: no grain formation (dotted 
line); grain formation with Gg eom — 7.07/iin (dashed curve); grain formation with a ge0 m = 0.707/um 
(solid curve). These three cases correspond respectively to models 1, 3, and 2 in Tables 1 and 2. 
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R = 0.5 AU, as well as at R = 1.0 AU (used in the previous figures). In all cases, the transitions 
occur in sequence at increasing vertical columns: first H to H2 at the low abundance of H2 (~ 10~ 6 ) 
needed to start the warm radical reaction sequence described in Sec. 2.1; then C to CO, next O to 
H2O, and finally full conversion of H to H2. Full conversion of the residual O to H2O occurs where 
x(H 2 ) ~0.1. 

We conclude that the formation of H2 on dust grains can play a key role in the synthesis of 
water in protoplanetary disk atmospheres. The efficiency of grain formation is determined overall 
by the effective surface area of the grains, which we have expressed in terms of the geometric 
mean a gC0 m of the minimum and maximum grain sizes in a MRN distribution. Both a ge0 m and the 
efficiency of grain formation are affected by the amount of grain growth and settling. Figure 2 gives 
a concrete illustration of this effect at 1 AU. The top panel of shows that increasing GWom from 
0.707/um to 7.07/im shifts the H-H2 transition downwards (to larger values of Nn). This means 
that the layer where the water abundance is a maximum is also shifted downwards, as shown in 
the dashed curves in the middle panel of Figure 2. 

We have focused on two aspects of grain formation of H2 that are important in disks, the 
extension of the standard theory to warm dust temperatures and the effect of grain growth and 
settling on the effective grain surface area. For the latter we use a single average parameter, the 
geometric mean of the minimum and maximum grain size (a geom ), ignoring many complexities such 
as the nature of the grain surface and the variation of dust properties with grain size. We also 
assume that the same grain area applies to both the formation of H2 and to the thermal coupling 
of the gas and the dust. The issue of H2 formation on warm grains arose earlier in connection with 
ISO measurements of the pure rotational transitions of H2 in photon dominated regions (PDRs) 
that seemed to require higher temperatures or more H2 than predicted by PDR models (e.g., van 
Dishoeck 2004). Habart et al. (2004) confronted this problem by considering the H2 grain formation 
rate as a variable in fitting measured 0-0 S(3) to 1-0 S(l) line ratios with a PDR model. They 
found "normal" (diffuse cloud) rate coefficients for high-excitation regions like the Orion bar and 
values enhanced by a factor of five for moderate excitation PDRs. They also determined a rough 
scaling between the H2 1-0 S(l) line and the strength of PAH features that suggested H2 formation 
on PAHS might be responsible for the empirically determined formation rates. Habart et al. (2004) 
and Cazaux et al. (2006) showed that "indirect" chemisorption, where a second H atom reaches 
a chemisorption site from a physisorption site, could explain the empirical grain formation rates 
deduced for PDRs if the barrier between the two kinds of sites is small enough (~ 0.05 eV) and if a 
significant population of small grains is present. While we do not consider H2 formation by PAHs 
here because there are likely to be few PAHs in the inner regions of protoplanetary disks (see the 
Appendix), we do employ a reduced population of small grains for the synthesis of H2. 

The question of whether a barrier as small a 0.05eV between physisorption and chemisorption 
sites is appropriate may have been resolved by laboratory experiments with graphite at moderate 
temperatures which suggest that the formation of H-atom dimers (Hornakasr et al. 2006a) or clusters 
(Hornakasr et al. 2006b) may be an important step in grain formation of H2. Although the Cazaux- 




Fig. 3. — : H2, CO and water abundances at R=0.25 (top), 0.5 (middle) and 1.0 (bottom) AU for 
model 2 with H2 formation on grains and a geom = 0.707/iin. 
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Tielens theory of grain formation of H2 yields a consistent picture for PDRs, it may not carry over 
to protoplanetary disk atmospheres of disks. The densities are much higher and the abundance 
of small grains is much less than in PDRs. However, the demonstrated sensitivity of molecular 
abundances to the abundance of H2 may eventually lead to a better understanding of H2 formation 
in disks when more extensive observations of species like water become available. 



2.5. The Formation of Water on Dust Grains 

Water may also be formed by the sticking of O atoms and OH radicals on warm grains and 
subsequent reactions with adsorbed hydrogen atoms or other species. In regions where water is 
not already the dominant oxygen species, atomic O is much more abundant than OH, and we 
consider only this case. For this purpose, we adapt the formula for H2 formation on grains (Eq. [To"]) 
by dropping the 1/2 factor, noticing that v(0) = (l/4)u(H) and introducing a grain formation 
efficiency factor e(O) and a sticking factor 5(0), 

i?(H 2 0) = 2 x 10- 19 cm 3 s- 1 T g 1 / 2 n(O)n H 

x(^)(^ M 0)S(0) (16) 
U.U1 a g 

We have implemented Eq. [16] using e(0) = 1 and S(0) = 1 to estimate the maximum role of water 
formation on grains. Figure 4 displays the changes in O and H2O for the case of no grain formation, 
i.e., pure gas phase synthesis (dotted curve), grain formation of only H2 with a ge0 m = 7.07/xm 
(dashed curve), and grain formation of both H2 and H2O with a ge0 m = 0.707^m (solid curve). 
These cases correspond respectively to models 1, 3, and 4 in Tables 1 and 2. This model of maximal 
grain formation leads to two effects for R = 1 AU. At very small columns, iVn < 2 x 10 20 cm -2 , the 
water abundance is increased by three dex to the 10 -10 level, but this level of hot water is probably 
not observable. Grain formation of water also plays some role deeper down (N-r > 10 21 cm -2 ) 
where gas phase synthesis is shut off by the low temperature. For Ah > 2 x 10 21 cm -2 at 1 AU, 
the increased abundance of water should hardly affect the amount of observable warm water since 
it occurs where the temperature is close to 100 K. However, inside of 0.5 AU, water at large depths 
Ah > 2 x 10 21 cm -2 will be observable because the gas is much warmer there. Nonetheless, the 
top layer of high- abundance water is the most important, and we conclude that grain formation of 
water is unlikely to play a dominant role in the inner disk. 

Another process that might contribute to the formation of water in the dense inner regions of 
a protoplanetary disk is radiative association, 

O + H 2 ^H 2 + hv. (17) 

To the best of our knowledge, this reaction has not been discussed in the literature. To exceed the 
production of water by grain formation, the rate coefficient for radiative association would have to 
be larger than 10 -17 cm 3 s _1 . Although this is a typical rate coefficient for the radiative association 



-15- 




Fig. 4. — : Oxygen species and temperature at 1 AU vs. vertical column density for three cases of 
H2 and H2O formation on grains: no grain formation (dotted line); H2 formation on grains with 
Ogeom = 0.707/im (dashed curve); - H2 and H2O formation on grains with cig eom = 0.707,™ (solid 
curve). These three cases correspond respectively to models 1, 3, and 4 in Tables 1 and 2. 
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of neutral species, it may not apply to water, at least without detailed consideration of the relevant 
potential energy curves. The rate coefficient for the radiative association of OH and H has been 
calculated by Field, Adams and Smith (1980). It is a strongly decreasing function of temperature, 
large at low temperatures (> 10 _14 cm 3 s _1 ) but small in warm regions (< 10~ 18 cm 3 s _1 above 
300 K). This process might be relevant for cool atomic H regions, but not in the atmospheres of 
the inner disks under consideration here. On the other hand, if the rate coefficient for Eq. [17] 
is significantly larger than 10~ 17 cm 3 s _1 at moderate temperatures, this reaction would be an 
important source of water. 



2.6. Mechanical Heating 

Finally we consider non-radiative heating processes that can affect the abundance and location 
of water. To illustrate this possibility, we adopt the mechanical heating model used in GNI04, 

9 

Tacc = |"h/5gC (18) 

where p g is the local gas mass density, c is the isothermal sound speed, Q is the angular rotation 
speed, and is a phenomenological constant. GNI04 used values in the range = 0.01 — 1.0; 
they showed that the value 0.01 differs negligibly from the case = for pure X-ray heating. 
Mechanical heating might arise from viscous accretion generated by the MRI or from the interaction 
of a stellar wind with the upper layers of the disk. The change from the dotted to the dashed curves 
in Figure 5 shows the effect of increasing from 0.01 to 1.0 for the case of no formation of H2 
on grains. As discussed in GNI04, mechanical heating causes the transition from high to low 
temperatures to start deeper down in the disk, near iVn = 10 21 cm -2 . It also produces a region of 
warm water (T > 300 K) at high abundance. The solid curve in Figure 5 shows that the combination 
of mechanical heating and grain formation of H2 can increase the thickness of warm water before 
gas phase production is cut off at low temperatures. The column of warm water for model 6 is now 
iV(H20) = 6.5 x 10 17 cm -2 at 1 AU, and even larger at smaller radii. These columns are of the 
same order as reported by Carr and Najita (2008) and Salyk et al. (2008). This result illustrates 
how mechanical heating can be important in determining the abundance of water and highlight the 
need for further studies of mechanical heating processes. 



3. Discussion 

In the previous section, we addressed the question of what determines observable levels of 
water in the inner region of a protoplanetary disk. One critical factor is the synthesis of H2, which 
is required for the production of OH and H2O from atomic O via neutral radical reactions. The 
gas must be warm in order for gas-phase synthesis to proceed, and the synthesis of H2 must occur 
at moderate depths into the transition region in order for the H + abundance to be reduced by the 
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1Q 20 1Q 21 1Q 22 1Q 

Vertical column density [cm -2 ] 



Fig. 5. — : Effect of mechanical heating on the abundance of molecular hydrogen and water: 
«h = 0.01 and no grain formation (dotted curve); = 1.0 and no grain formation (dashed 
curve); = 1.0 with grain formation of H2 and 

^geom — 0.707/xm (solid, curve). These three cases 
correspond respectively to models 1, 5, and 6 in Tables 1 and 2. 
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attenuation of the X-rays and by the increase in density. All of these factors, the H2 abundance, 
the temperature, and the level of ionization, can come together to produce maximum abundances 
of water in the thermal-chemical transition region that lies between the hot upper layer and the 
mid-plane of the atmosphere of an X-ray irradiated protoplanetary disk. 

We arrived at this conclusion by showing how sensitive the amount of warm water in proto- 
planetary disk atmospheres is to poorly known processes involving the H2 molecule. We used the 
device of comparing the H2 and H2O abundances in the original model of GNI04 (corrected for the 
error in the rate coefficient k§ in Eq. [5|) with the abundances calculated for new rate coefficients and 
processes. These include (1) the adoption of the Flower and Harris (2007) rate coefficients for the 
three-body reactions Eq. [31 (2) the introduction of grain formation of H2 on the basis of the theory 
of Cazaux and Tielens (2003, 2004), Eq. I15|. (3) the assumption that H2O is formed on grains at the 
maximum rate, Eq. [TBI and (4) mechanical heating of the atmosphere, Eq. [T5J These changes were 
illustrated by the abundance patterns in the previous figures. Although we emphasized water, the 
changes we have discussed also affect atomic O and OH and the species with which they interact. 




1Q 20 1Q 21 1Q 22 1Q 

Vertical column density [cm -2 ] 



Fig. 6. — : Water abundance and temperature profiles at R=0.25 (dashed), 0.5 (dotted) and 1.0 
(solid) AU for model 1 for the case of only gas phase reactions (model 1) and a geom = 0.707/um. 

An important measure of the changes is the column density of warm water, defined to have 
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T > 300 K and an abundance large enough to potentially make a significant contribution to 
observable line fluxes. In practice, this column begins where the water abundance first achieves 
its maximum value and ends where the temperature drops below about 300 K. These columns are 
given in Table 2 for the models considered in the previous section, along with the columns of 
warm OH and CO. The table includes the corresponding ranges in hydrogen column density (AT H ), 
temperature (T), electron fraction (x c ), and molecular hydrogen abundance (x(H2)). Not only 
do these abundances express the physical conditions where the warm water occurs, but they also 
determine which collision partners might be important in exciting observable transitions. Electronic 
excitation of the low-lying rotational and ro- vibrational levels of water by electrons are typically 10 4 
times stronger than for H2 (Faure and Josselin, 2008). Thus electronic excitation will not dominate 
in these warm water regions. But the H-H2 transition may not be complete, so that atomic H 
collisional excitation is likely to play a role, not just H2 collisions. 

Table 2 contains information about the radial variation of the warm molecular columns of 
H2O, OH and CO. The warm water column iV(H20) rapidly increases with decreasing radius, 
faster than 1/R for R < 0.5 AU. The inner disk gets warmer due to increased accretion heating, 
which reduces the cooling of the gas by the dust. All of the models yield iV(H20) > 10 18 cm~ 2 for 
small radii (R < 0.5 AU). This result opens up the possibility that the water lines may have broad 
wings for inclined disks. Detailed excitation calculations are needed to deal with this issue, not 
just for water but for other molecules for which line shape measurements are feasible, especially 
CO. Accretion heating becomes increasingly important with decreasing R, to the point where pure 
gas phase processes (model 5 in Table 1) give about the same water columns as model 6 (with H2 
grain formation). This can be seen in detail in Figure 6 which shows water abundance profiles at 
several radii for pure gas phase chemistry. 

The results of the few reported observations of water in disks around young low-mass stars 
(Carr et al. 2004; Carr and Najita 2008; Salyk et al. (2008) have been expressed in terms of a column 
density and an excitation temperature. For the typical T Tauri star AA Tau, Carr and Najita find 
iV(H20) = 6.5 x 10 17 cm~ 2 and a characteristic excitation temperature of ~ 575 K. The same 
quantities for OH and CO are iV(OH) = 8.1 x 10 16 cm~ 2 (characteristic temperature of ~ 525 K) 
and 7V(CO) = 4.9 x 10 17 cnr 2 (characteristic temperature of ~ 900 K). Salyk et al. (2008) report 
similar water columns for two other T Tauri stars. These column densities have been derived 
assuming thermal populations, so they may be underestimates or overestimates, depending on 
the excitation conditions. Although the model calculations in this paper were done primarily to 
illustrate how certain physical and chemical processes affect the amount of observable water in 
protoplanetary disks - and not to model any particular disk, it is of some interest to compare 
the observations with the results in Table 2. The presumed improvements in the treatment of the 
formation of H2, especially grain formation, bring the calculated column of warm water close to 
the observed range iV(H20) > 10 17 — 10 18 cm~ 2 . If the heating of the atmosphere is increased, 
the reported columns iV(H20) = 6 — 8 x 10 17 cm~ 2 can be reached or even exceeded. In drawing 
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Table 2:: Warm molecular columns 



R=0.25 AU 


Model 


JV H [cm -2 ] 


T[K] 


x(e) 


x(R 2 ) 


iV(H 2 


) [cm- 2 ] 


JV(CO) [cm- 2 ] 


N (OH 


[cm" 2 ] 


1° 


1.7(21)-9.5(21) 


300-560 


2.0(-6)-4.2(-6) 


5.5(-2)-4.5(-l) 


1.1 


18) 


1.1(18) 


2.8 


15) 


2 b 


9.3(20)-9.5(21) 


300-810 


1.6(-6)-7.9(-6) 


7.3(-2)-4.9(-l) 


1.5 


18) 


1.1(18) 


2.4 


15) 


3 C 


1.7(21)-1.4(22) 


310-650 


1.5(-6)-4.2(-6) 


4.1(-2)-4.9(-l) 


2.8 


18) 


1.6(18) 


1.3 


15) 


4 d 


9.3(20)-9.5(21) 


300-810 


1.5(-6)-7.9(-6) 


7.5(-2)-4.9(-l) 


1.6 


18) 


1.1(18) 


2.4 


15) 


5 e 


2.1(21)-1.9(22) 


300-720 


7.4(-7)-4.1(-6) 


4.7(-2)-4.9(-l) 


3.7 


18) 


2.3(18) 


1.4 


15) 


6' 


1.4(21)-1.9(22) 


300-1000 


5.7(-7)-l.l(-5) 


3.4(-l)-4.9(-l) 


4.0 


18) 


2.3(18) 


7.9 


14) 


R=0.5 AU 


Model 


JVh [cm -2 ] 


T[K] 


x(e) 


x(K 2 ) 


iV(H 2 


) [cm- 2 ] 


AT(CO) [cm- 2 ] 


N (OH 


[cm" 2 ] 


1° 


1.2(21)-2.1(21) 


320-510 


4.7(-6)-5.7(-6) 


4.4(-2)-2.6(-l) 


6.2 


15) 


1.0(17) 


7.3 


14) 


2 b 


4.3(20)-2.1(21) 


320-790 


5.1(-6)-9.2(-6) 


2.9(-2)-4.7(-l) 


3.0 


17) 


2.0(17) 


1.7 


15) 


3 C 


9.2(20)-4.8(21) 


330-600 


2.0(-6)-3.5(-6) 


3.5(-2)-4.7(-l) 


7.9 


17) 


4.9(17) 


1.0 


15) 


4 d 


5.5(20)-2.1(21) 


320-790 


5.0(-6)-8.5(-6) 


3.1(-2)-4.7(-l) 


3.1 


17) 


2.0(17) 


1.7 


15) 


5 e 


1.9(21)-6.4(21) 


300-670 


1.6(-6)-2.8(-6) 


6.6(-2)-4.4(-l) 


9.2 


17) 


6.5(17) 


8.8 


14) 


6' 


6.2(20)-5.5(21) 


320-1090 


1.5(-6)-l.l(-5) 


2.0(-l)-5.0(-l) 


1.1 


18) 


5.6(17) 


7.9 


14) 


R=l AU 


Model 


JV H [cm" 2 ] 


T[K] 


x(e) 


x(H a ) 


iV(H 2 


) [cm- 2 ] 


JV(CO) [cm- 2 ] 


N (OH 


[cm" 2 ] 


1° 


> 1.2(21) 


>300 


>1.2(-5) 


<7.5(-3) 


4.4 


12) 


6.6(15) 


6.2 


13) 


2" 


4.3(20)-1.2(21) 


330-770 


6.3(-6)-8.1(-6) 


9.7(-2)-4.6(-l) 


1.6 


17) 


9.9(16) 


7.9 


14) 


3 C 


1.0(21)-2.9(21) 


310-560 


2.7(-6)-3.3(-6) 


6.4(-2)-4.5(-l) 


4.0 


17) 


2.9(17) 


7.9 


14) 


4 d 


4.3(20)-1.2(21) 


330-770 


6.2(-6)-8.1(-6) 


9.9(-2)-4.6(-l) 


1.6 


17) 


9.2(16) 


6.6 


14) 


5 e 


1.8(21)-3.8(21) 


340-630 


1.6(-6)-2.1(-6) 


7.4(-2)-3.7(-l) 


4.5 


17) 


3.7(17) 


5.1 


14) 


6' 


1.0(21)-3.8(21) 


300-1020 


1.6(-6)-8.0(-6) 


4.6(-l)-5.0(-l) 


6.5 


17) 


3.3(17) 


1.6 


14) 


R=2 AU 


Model 


JVh [cm -2 ] 


T[K] 


x(e) 


x(H a ) 


iV(H 2 


) [cm- 2 ] 


JV(CO) [cm- 2 ] 


N (OH 


[cm" 2 ] 


1° 


> 9.4(20) 


>300 


>3.7(-5) 


<l.l(-3) 


3.0 


11) 


5.0(13) 


3.2 


13) 


2 ft 


3.9(20)-9.4(20) 


300-720 


6.6(-6)-8.6(-6) 


1.4(-l)-4.4(-l) 


1.0 


17) 


8.4(16) 


1.1 


15) 


3 C 


9.4(20)-1.8(21) 


310-500 


3.4(-6)-3.7(-6) 


1.0(-1)-4.1(-1) 


1.7 


17) 


1.7(17) 


6.2 


14) 


4 d 


3.9(20)-9.4(20) 


300-720 


6.6(-6)-8.6(-6) 


1.4(-l)-4.5(-l) 


1.0 


17) 


8.4(16) 


1.2 


15) 


5 e 


1.8(21)-3.4(21) 


320-600 


1.4(-6)-1.9(-6) 


9.8(-3)-3.4(-l) 


3.0 


17) 


3.2(16) 


6.2 


14) 


6' 


4.5(20)-2.9(21) 


310-1170 


1.7(-6)-l.l(-5) 


3.0(-l)-5.0(-l) 


5.7 


17) 


2.8(17) 


4.4 


14) 


a Gas ph 


ase only, a geom = 


0.707/im, a 


= 0.01 















6 H 2 formation on grains, a geom = 0.707/im, a = 0.01 

C H 2 formation on grains, a georn — 7.07/im, a = 0.01 

d H 2 and H 2 formation on grains, a geom = 0.707/im, a = 0.01 

e Gas phase only, a geom = 0.707/im, a = 1.0 

■^H 2 formation on grains, a geom = 0.707/im, a = 1.0 
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this conclusion, we have not varied the density distribution or the X-ray irradiation, which can 
affect the water column. For example, the X-ray emission of AA Tau appears to be typical for an 
active classical T Tauri star (Grosso et al. 2007), with a minimum or quiescent X-ray luminosity 
Lx = 5 x 10 30 ergs _1 and temperature Tx = 2.5 keV. Both parameters are larger than used here 
based on the GNI model, Lx = 2 x 10 30 ergs _1 and temperature Tx = 1.0 keV. In addition to the 
effects of choices in model parameters, the predictions in Table 2 are subject to the limitations of 
our demonstration model (GNI04), e.g., the deficiencies in the thermal model and specifically the 
neglect of water heating and cooling, as discussed in Sec. 2.1. 

One feature of the model calculations is that the column density of OH is much less than that 
of water. According to the last column of Table 2, the typical column density of warm OH in the 
inner disk is, to within a factor of two, A r (OH) ~ 10 15 cm -2 . The ratio of the warm H2O to the 
warm OH column density ranges from a few hundred to several thousand, due mainly to the changes 
in the chemical model. The measured warm OH column for AA Tau is iV(OH) = 8 x 10 16 cm -2 , 
and none of the entries in Table 2 approach this value. This points to a deficiency in our calculation 
of the OH abundance. Staying within the confines of gas-phase chemistry, this problem might be 
solved by including FUV photo-dissociation of OH and H2O, taking into account that the cross 
section for water is several times larger than for OH. Preferential photo-dissociation of water might 
then lead to a significant increase in the amount of OH. An extreme example of this effect is the 
detection in the L band (3/im) of the lines of OH but not H2O in the disks around two Herbig Ae 
stars (Mandell et al. 2008). In this case there is ample FUV from the A star to destroy circumstellar 
water in favor of OH. For T Tauri star disks, however, the stellar FUV flux is much smaller, and it 
may be unable to shift the balance from H2O to OH, especially in the presence of neutral gas phase 
synthesis reactions which tend to drive OH to H2O. This is the situation in our model, despite the 
fact that the rate coefficient for the reaction that destroys water, Eq. [9] (H + + H2O —* H20 + + H; 
Anicich 1993), is four times larger than that for OH, Eq. M (H + + OH -> OH+ + +H; UMIST guess, 
Woodall et al. 2007). However, Bergin et al. (2003) have pointed out that the photo-dissociation 
of water is sensitive to the spectral shape of the stellar FUV, in particular to the strength of the 
Lyman-a line, which dominates the FUV emission in some T Tauri stars such as TW Hya. Thus 
sources with stronger Ly-a emission would be more effective in destroying water. On the other 
hand, OH can also be destroyed by Ly-a radiation (van Dishoeck and Dalgarno 1983). The rates of 
direct destruction of OH and H2O by X-rays are about the same, since both cases depend largely 
on K-shell absorption by a single oxygen atom. 

Going beyond pure gas-phase chemistry, OH may be formed by desorption from icy grain 
mantles transported from cold regions to the warm source of the observed MIR lines within an 
AU or so from the star. The release could proceed by either thermal or photo-desorption (Ciesla 
and Cuzzi 2006). The latter process has been suggested to be the origin of the highly excited OH 
ground rotational lines seen in HH 211 (Tappe et al. 2008). In this case, the FUV radiation is 
shock generated, and it can desorb both OH and H2O (e.g., Andersson and van Dishoeck 2008). 
When the water is photo-dissociated, OH is formed in a highly excited state, and the excited OH 
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rotational lines are produced by fluorescence. Similar processes may be at work in T Tauri disks, 
mediated by stellar FUV or X-ray radiation. 

A last comparison between the present model calculations and the limited observations to test 
them involves CO. According to Table 2, the model calculations predict comparable warm columns 
of H2O and CO. The average ratio of water to CO to water columns in Table 2 is 1.6, independent 
of radius, with small variations of order 25%. Carr and Najita (2008) report a ratio of 3/4 for 
AA Tau, in rough accord with our model which ascribes most of the warm water and CO to the 
thermal-chemical transition region. The higher temperature and smaller emitting area of the CO 
lines compared to the water lines in AA Tau suggest that the CO and the water probe somewhat 
different radial distances. While the CO likely arises from small radii (R < 1 AU), the water may 
arise from radii out to > 2 AU (Carr & Najita 2008). The individual column densities of warm CO 
and water reported in Table 2 are consistent with the properties of the molecular emission reported 
for AA Tau. 

It is also of interest to compare our results with other chemical calculations of the inner disk, 
e.g., Markwick et al. (2002), Agundez et al. (2008), Woods and Willacy (2008), and Gorti and 
Hollenbach (2008). These comparisons are characterized by the fact that disk chemistry models 
differ widely in their underlying assumptions as well as in execution. None of these studies include 
mechanical heating, and only Woods and Willacy consider the recent theory of grain formation 
of H2 by Cazaux and Tielens (as do Nomura et al. 2009). For example, Markwick et al. (2002) 
report only total molecular column densities with no information about the vertical abundance 
distribution. At 1 AU, they obtain N(H 2 0) = 1.6 x 10 22 cm -2 and 7V H = 1.6 x 10 26 cm' 2 , for 
an average water abundance x(H20) = 10~ 4 . Since they include only viscous heating, their disk 
atmosphere is much cooler than obtained in models that include stellar irradiation, and most of 
their calculated water column will not be warm. Markwick et al. do not discuss OH. 

Agundez et al. (2008) developed a time-dependent model oriented towards the synthesis of 
complex molecules with neutral radical reactions. In Sec. 4 of their paper, they use the D'Alessio 
(2001) model exposed to stellar and interstellar FUV radiation, but without X-rays. They assume 
equal gas and dust temperatures and start with a molecular mixture with all oxygen in water that is 
not in CO. They obtain the following steady-state columns inside 3 AU, approximately independent 
of radius: iV(CO) = 10 18 cm" 2 , N(R 2 0) » 10 17 cm" 2 and iV(OH) = 10 14 cm" 2 . Although these 
values are similar to those calculated here, the molecules will not be as warm nor produce as much 
emission because the gas temperature is the same as the dust temperature. 

Woods & Willacy (2008) have developed a thermal-chemical model based on the D'Alessio 
(2001) disk structure that includes stellar and interstellar FUV radiation and X-rays. Their main 
emphasis is on the carbon isotopes, but they also report on the total water column density and give 
contours of water and OH in the inner disk. They show that the regions of the disk atmosphere, 
where the average temperatures of the OH and water emission are the same as found for these 
species in AA Tau (Carr and Najita 2008), have local abundances in the model that are similar to 
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those reported by Carr and Najita. No quantitative information is provided on the column densities 
of warm OH and water. 

Gorti and Hollenbach (2008) have calculated a wide variety of atomic and molecular lines 
with a full thermal-chemical model that is similar to but more general than the one used here. It 
includes X-rays as well as stellar FUV, although the X-ray heating is significantly less than in our 
model. Gorti and Hollenbach also include cosmic rays, whereas we believe, by analogy with the 
much less powerful solar wind, that they are blown away from the inner protoplanetary disk by 
the stellar wind of the YSO (Glassgold, Najita and Igea, 1997). A key element of their model is 
its use of a significant population of small grains and PAHS to heat the gas via the photoelectric 
effect. However, they do not include mechanical heating, which we find plays an important role 
in the inner disk. Gorti and Hollenbach also simultaneously solve the equation for hydrostatic 
equilibrium, so that their disk atmosphere is more flared than ours. The column density of water 
is smaller than calculated here, and the line strengths are modest because the warm water is at 
smaller densities due to disk flaring (D. Hollenbach, private communication 2009). Consequently, 
they predict that only a few MIR lines should be detectable by the Spitzer IRS, in contrast with 
the rich emission spectra reported by Carr and Najita (2008) and Salyk et al. (2008). On the other 
hand, we do not treat hydrostatic equilibrium here, and our future modeling will need to include it 
in order to assess its importance for the formation of H2 and the generation of significant columns 
of warm water. 



4. Conclusions 

Our results indicate that the conditions that lead to strong water emission in the infrared, i.e., 
large column densities of warm water molecules near the disk surface, can be produced in situ under 
certain conditions in an X-ray irradiated disk model without requiring radial or vertical transport 
of water to the disk surface. In the inner disk (< 0.25 AU), abundant warm water can be produced 
through gas-phase processes and only X-ray heating. Efficient H2 formation on warm grains and 
non-radiative heating appear to be capable of enlarging the radial range over which warm water 
emission is produced. These processes, when taken together, can account roughly for the column 
densities and emitting areas of the MIR water emission reported for T Tauri disks by Carr & Najita 
(2008) and Salyk et al. (2008). 

We find that the column density (and temperature) of water at the disk surface depends 
significantly on the grain surface area in the disk atmosphere (i.e., on the extent of grain growth and 
settling) and on the heating of the atmosphere (and therefore perhaps on the disk accretion rate). 
Because both the disk accretion rate and the extent of grain settling are believed to vary among 
T Tauri disks, our model results suggest that the strength and character of MIR water emission 
from disks may also show some diversity. In contrast, because much less restrictive conditions are 
needed to produce large column densities of warm water close in (< 0.25 AU), water emission from 
this inner region of the disk may be quite common. 
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More specifically, for the grain areas that we use, we find that sources with modest amounts of 
grain growth (e.g., a factor of 20 decrease in area compared to the ISM) can produce large column 
densities of warm water within 1 AU. Even larger columns result for ten times more grain growth 
(see Table 2), and we might expect this trend to continue for a growth factor of 1000. Recent studies 
of the MIR spectral energy distributions of T Tauri stars suggest that disks have experienced grain 
growth factors spanning a similar range (factors 10-1000; Furlan et al. 2006) as those considered in 
our models. Thus, we might expect that T Tauri disks that show different amounts of grain growth 
may display a diversity in the strength and character of their water emission. 

Similarly, T Tauri disks are also believed to have accretion rates that vary by at least two orders 
of magnitude (Hartmann et al. 2006). Stellar X-ray luminosities also vary significantly. Since 
these factors influence the ionization and amount of surface heating, they may also contribute 
additional diversity to the water emission spectra of T Tauri disks. Large accretion rates can 
enhance the heating in the atmosphere and increase the strength of the water emission. A higher 
X-ray luminosity would produce both stronger ionization and heating. According to the simplified 
theory of Sec. 2.2, heating effects are more important than the ionization level. Such trends may be 
revealed in Spitzer IRS data sets that are being collected currently. However, if radial or vertical 
mixing is efficient, it might diminish or erase the correlations that would be predicted by our models. 

In addition to searching for correlations in large data sets, it may also be interesting to explore 
individual objects in greater detail in the context of our models. In the case of AA Tau, the analysis 
of Furlan et al. (2006) suggests that the dust in the disk surface may be depleted by a factor of 
nearly 1000. Our calculations suggest that the observed warm water columns can be accounted 
for, even for this degree of depletion without non-radiative heating. Thus in this case, mixing and 
transport appear to be unnecessary to account for the observations. 

We have identified two major uncertainties in the chemistry of H2 that bear on the amount of 
water and other species: (a) Three-body formation 0/H2 - Although the use by Flower and Harris of 
a new equilibrium constant represents an important step forward, the data on which these reactions 
are based are old measurements of the low-energy collisional excitation of H2 by atomic hydrogen 
(Jacobs et al. 1967; Baulch et al. 2005). New experiments are needed, (b) Grain formation 0/H2 - 
The adopted rate coefficients are based on theory (Cazaux and Tielens 2002, 2004). Although this 
model is consistent with observations of H2 lines in PDRs (Habart et al. 2004, Cazaux and Tielens 
2006), its applicability to the dense and dust-depleted atmospheres of protoplanetary disks is not 
yet fully established. 

Regarding non-radiative heating, we have used Eq. [18] from GNI04, who suggested two ways 
of heating the transition region above the level achievable with X-rays, dissipation of turbulence 
and the interaction of the stellar wind with the disk atmosphere. A sound basis for these processes 
has not yet been laid, but progress should be possible with improved numerical simulations of the 
magneto-rotational instability (e.g., Turner and Sano 2008; Terquem 2008) and of the wind-disk 
interaction (e.g., Li and Shu 1996; Matsuyama et al. 2009). Our results highlight the need for 
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further basic studies of the heating generated by the MRI and by the interaction of the wind and 
the disk. 

Given these current uncertainties, and the potential sensitivity of water emission to multiple 
effects (X-ray irradiation, grain growth, disk accretion, transport, etc.), it may be difficult to 
diagnose the role of such important (but elusive) processes as radial and/or vertical transport from 
observations of water alone. It would be worthwhile to establish diagnostics for other relevant 
regions of a protoplanetary disk such as the disk midplane and the outer disk. An interesting 
example is provided by the discussion by Woods and Willacy (2007) of the formation of benzene 
near the midplane of the inner part of a protoplanetary disk. 
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A. PAH Formation of H2 in Protoplanetary Disks 

Duley and Williams (1993) suggested that the H2 forms on PAHs in carbonaceous interstellar 
dust, and Bauschlicher (1998) demonstrated that it was energetically favorable for PAH cations 
to acquire an H atom and then form H2 in a subsequent collision with atomic H. He suggested 
that this was a way to synthesize H2 in the gas phase. Rauls and Hornekaer (2008) have recently 
discussed how super-hydrogenated neutral PAHs might also play a role in the gas-phase synthesis 
of H2. Jonkheid et al. (2006) included H2 formation by PAHs in modeling the transitional disk 
around HD 141569. Habart et al. (2004) adduced possible evidence for H2 formation by PAHs in 
PDRs. 

As a basis for evaluating the possibility that PAH formation of H2 is relevant in the molecular 
transition region of disks, we make a crude estimate of the rate coefficient for the Bauschlicher 
process, 

(PAH + + H) + H^PAH+ + H 2 , (Al) 

where the symbol (PAH + + H) represents a stable PAH cation with an added H atom. The effective 
rate coefficient -Rpah for this process, which can then be compared with the rate coefficient R for 
grain formation in Eq. [T5l can be written as, 

Rpar = A:pah /(PAH+ + H) xpah, (A2) 

where /cpah is the unknown rate coefficient for reaction IA1|. /(PAH + + H) is the fraction of 
appropriately hydrogenated PAH cations, and xpah is the total gas phase abundance of PAHs 
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in all ionization states. We make a rough guess that /cpah ~ 10~ 10 cm 3 s _1 and assume that 
/(PAH + + H) ~ 10 -2 and xpah ~ 10 -7 , an upper limit to the range deduced by Geers et al. (2006) 
from Spitzer observations of disks. The result, -Rpah ~ 10~ 19 cm 3 s , is only slightly smaller than 
given by Eq. [15] for grain formation of H2 . 

Each of the factors in our estimate of -Rpah is uncertain by at least one dex, so we cannot rule 
PAH formation of H2 in or out on this basis. However, we suspect that it is unimportant for the 
inner regions of protoplanetary disks because PAHS close to YSOs are likely to be destroyed by 
stellar X-rays. We estimated the total abundance xpah ~ 10~ 7 on the basis of the low detection 
frequency of PAH IR features in Spitzer observations of T Tauri stars (Geers et al. 2006; ~ 10 %) 
and Class I YSOs (Geers et al. 2008; < 2%). Even where observed, the PAHs may occur mainly at 
large radii, as for Herbig stars (Habart et al. 2006; Geers et al. 2007; Goto et al. 2008). Laboratory 
experiments with the Synchrotron Radiation Facility in Grenoble show that PAHs and very small 
particles are destroyed by hard X-rays (Guegeon 1998; Mitchell et al. 2002; see also Sec. 6 of the 
review by Glassgold, Feigelson & Montmerle 2000). Indeed, Voit (1992) invoked X-ray destruction 
of PAHs to explain the difference between the observed PAH features seen in starburst galaxies and 
AGN. Finally, Siebenmorgen and Kreugel (2008) have recently modeled the PAH abundance in T 
Tauri disks and find that energetic photons destroy PAHS inside of 10 AU. It is for these reasons 
in support of a low PAH abundance that we tentatively ignore PAH formation of H2 for the inner 
part of the atmospheres of protoplanetary disks. 
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